Anomalous Hall Effect due to the spin chirality in the Kagome lattice 
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We consider a model for a two dimensional electron gas moving on a kagome lattice and locally 
coupled to a chiral magnetic texture. We show that the transverse conductivity a xy does not vanish 
even if spin-orbit coupling is not present and it may exhibit unusual behavior. Model parameters 
are the chirality, the number of conduction electrons and the amplitude of the local coupling. Upon 
varying these parameters, a topological transition characterized by change of the band Chern num- 
bers occur. As a consequence, a xy can be quantized, proportional to the chirality or have a non 
monotonic behavior upon varying these parameters. 
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INTRODUCTION 



In ferromagnetic systems, there are two contributions 
to the transverse resistivity p xy : one is due to the usual 
Lorentz force acting on the electrons when a magnetic 
field is applied, RqB, the second one, R S M, is propor- 
tional to the magnetization of the ferromagnet. This is 
called the anomalous Hall effect. 

The origin of this anomalous Hall effect has long been 
controversial and both extrinsic (impurities) or intrin- 
sic mechanisms were discussed. Karplus and Luttingesi 
proposed that this effect is a consequence of spin-orbit 
interaction in metallic ferromagnets. Then it was argued 
that impurities give the main contribution to the anoma- 
lous Hall effect and usually two mechanisms, both due 
to spin-orbit coupling, contribute: one is known as side- 
jump mechanism 2 *^ and it predicts that the Hall resis- 
tivity R s is proportional to the longitudinal resistivity 
p. The second one 4 ^ is the skew-scattering mechanism, 
which gives a contribution proportional to p 2 . 

In the recent years several groups measured anomalous 
Hall effect in various systems which cannot be attributed 
to usual mechanisms (skew-scattering or side-jump). A 
new intrinsic mechanism, related to non-collinear spin 
configuration, with a ferromagnetic component, was first 
proposed for manganite ofi'Tifi : in these systems deviation 
from collinearity is a consequence of the competition be- 
tween double exchange, superexchange and spin-orbit in- 
teractions. Then it was proposed that a similar mecha- 
nism works in spin glass systems where spin configuration 
is highly non coplanauSiifl: in the weak coupling limit it 
was shown that the Hall resistivity is proportional to the 
uniform chirality parameter with a sign which depends on 
the details of the band structure. This chirality parame- 
ter was defined by Tatara and Kawamurai for three spins 
at sites i, j, k as a scalar quantity: \ — (S».(Sj x Sfe)), 
being non-zero if the spins are non-coplanar. Very re- 
cently several groups confirmed the existence of such a 
contribution in typical AuFe spin glasses!!* 1 ^. 

In non-disordered systems, the same mechanism can 
work as soon as the ordered magnetic structure is non- 
coplanar. This is realized in several pyrochlore com- 
pounds such as Nd2Mo2C>7 for which Taguchi et alii pro- 



posed that an anomalous contribution to the Hall effect is 
related to the umbrella structure of both Nd and Mo mo- 
ments in the long range ordered phase. Detailed studies 
of this system and other pyrochlores (Sm 3; Yi_ a ;Mo207, 
Nd x Yi_ 2 ;Mo207) have been performe d 1 ^ 14 ' 1,1 ? leading to 
some controversy: using neutron diffraction experiments 
in magnetic field it is possible to calculate the T and H 
dependence of both Mo and Nd chiral order parameter 
and in Ref. it was concluded that the chiral mecha- 
nism alone cannot explain the T-dependence of the Hall 
coefficient in pyrochlores. 

From the theoretical point of view, this mechanism is 
often called "Berry phase contribution" because a non- 
vanishing spin chirality is associated with a non- vanishing 
Berry phase for conduction electrons coupled via local 
Hund's exchange interaction to the spins. Ohgushi et 
alii calculated the Hall effect in a Kagome lattice with a 
non-coplanar long range spin structure. They have shown 
that, in the adiabatic limit, when the electron conduc- 
tion spins are aligned with localized spins at each site 
of the lattice (this corresponds to infinite Hund's cou- 
pling) , the Berry phase contribution to Hall conductivity 
is quantized for some values of the band filling. In this 
paper, we study a model which extrapolates between this 
strong coupling limit and the weak coupling case studied 
by Tatara and Kawamura®. We show also that the Berry 
phase contribution does not depend only on the chirality, 
but also on the strength of the local Hund's coupling and 
on the band filling. 



II. MODEL 

A. Hamiltonian 

The aim of this work is to show how the transport prop- 
erties of electrons are influenced by two different contri- 
butions. First, the electrons are restricted to move on a 
lattice and are therefore experiencing its geometry. This 
results in a peculiar band structure. Second, each elec- 
tron has its spin locally coupled to a given distribution 
of magnetic moments on each site of the lattice. 
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Both effects are taken into account in the following 
Hamiltonian: 

U = { c i* c i° + h - c ) ~ J Y c \ a {°<*p ' Si) (Hp. 

(1) 

The first term describes the electrons moving on the lat- 
tice: tij is the hopping integral between two neighboring 
sites i and j; c\ a and c ia are the creation and annihilation 
operators of an electron with spin a on the site i. The 
second part of the Hamiltonian couples the electron spin 
to a local moment on each site. 

The coupling constant to each local moment Si is J, 
and these moments are treated below as classical vari- 
ables. (Tap are the Pauli matrices. The underlying lat- 
tice is the Kagome lattice, depicted in Fig.Q] It is a two 
dimensional tiling of corner sharing triangles, described 
as a triangular lattice of triangles throughout this article. 

This Hamiltonian has already been discussed in the 
limit of J — > oo by Ohgushi et alil& In this limit the 
two a =t,4 bands are infinitely splitted and the model 
describes a fully polarized electron gas subject to a mod- 
ulation of a fictitious magnetic field, corresponding to 
the molecular field associated with the magnetic texture. 
The comparison with the present work will be done later 
on. In a general case of finite J, the calculations must be 
done numerically. It is worth noting that the sign of J is 
unimportant in this classical treatment since changing J 
to — J is equivalent to an exchange of \ and J. spin states. 




FIG. 1: (a) The Kagome lattice is described as a triangu- 
lar lattice of triangles. Its Bravais vectors are a\ = (1,0), 
a-z = ( — 1/2, \/3/2), and a- A — —(1/2, \/3/2), connecting re- 
spectively, the points A to B, B to C, and C to A. (b) The 
first Brillouin zone is an hexagon with the corners located at 
k = ±(2tt/3) ai, k = ±(2tt/3) a 2 and k = ±(2tt/3) as- (c) The 
umbrella structure on the triangular cell of the kagome lattice. 



we are interested in describing the anomalous trans- 
port properties of the t ype already addressed in previ- 
ous works (see Refs. [l^llflllTl p in which it has been ar- 
gued that the transverse conductivity may be related to 
the chirality of the magnetic phase. Second, an unusual 
mechanism of the AHE has recently been proposed, based 
on an unconventional chiral magnetic ordering due to ge- 
ometrical frustrationi^iSiiSiSfi but with controversial in- 
terpretations. We therefore consider the present model 
as a simple model to address these problems and clearly 
identify each relevant contribution. 

We consider the magnetic phase which has been 
obtained by studying a pure spin model with the 
anisotropic Dzyaloshinskii-Moriya interactions on the 
Kagome lattice^!. It consists of an umbrella of three spins 
per unit cell of the Kagome lattice (see Fig.^J. Each um- 
brella can be described by the spherical coordinates of the 
three spins (n/6, 9), (5tt/6, 9) and (— 7r/2, 9). The angle 
9 ranges from (all the spins are perpendicular to plane, 
and the corresponding ordering is ferromagnetic) to n. 
It is worth noting that the usual three sublattice planar 
magnetic phase belonging to the ground state manifold 
of the kagome antiferromagnet, and commonly denoted 
as the "q — phase" corresponds to 9 = ir/2. 

In between, all phases are possible and chiral, where 
chirality is defined as the mixed product of three spins 
on a plaquette 

Xijk = St ■ (S, x S fc ) = ^ cos 9 sin 2 9. (2) 

All these phases (9 e [0, n]) are therefore translation- 
invariant but do not have time reversal symmetry. 

We consider the local spins as classical spins with 
length M and we introduce in eq.QJan effective coupling 
constant Jo = JM which allows to rewrite the hamilto- 
nian as 

U = ( c i<T C 3>? + h - c ) ~ J o C L {Veep ■ n») (kp 

{i,i),<r « 

(3) 

where is a unit vector collinear to the local moment 
Si. 

Setting t = \tij\ as the energy unit, our free param- 
eters are the angle 9, parametrizing the chirality of the 
magnetic texture, and the value of Jo . The Fermi level, 
through the filling factor p of the bands, can also be var- 
ied. Each of these variables gives a different contribution 
as will be discussed in the next sections. 



B. Parameters and motivation 

We consider the spin texture as a 'parameter' of the 
model. The local moments Si are described classically, 
and the chosen magnetic phase is periodic, thus allowing 
to work in the reciprocal space. The choice of the mag- 
netic arrangement is motivated by two reasons. First, 



C. Physical quantities 

Before calculating any physical observables, we have to 
diagonalize the Hamiltonian of Eq. For this purpose, 
the Hamiltonian is rewritten in the reciprocal space as 

W = Z)*kfck*k + h.c, (4) 

k 
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with \l/k = (cAkT, CBkT, cckt, c A kU c Bki, Cckl), A -B and 
C are the corners of the kagome lattice unit cell (Fig. ^ 
(a)) and 

CAko = c A,j,« e ik ' T ' 3 ', (5) 
j 

I 



is the Fourier transform of CA.j,cr and r^j = ta + -Rj with 
i?j the lattice vector and ta the position of the moment 
Sa in the unit cell, hk is a 6 x 6 matrix given by 



hk = 
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(6) 



with ty = i and = 2vicos(k • a,). 

To calculate the Anomalous Hall Effect, we use the 
expression of the off-diagonal conductivity using the fol- 
lowing Kubo formula in the limit of disorder-free electron 
gas2£ 



VxyW) - — ^ ^ (fnk " fmk) 
n^m k 

(n, k\v x \m, k) (n, k\v y \m, k) 



(£nk — £mk) (^nk ~ £mk — w) ' 



(7) 



where S is the surface of the unit cell, e n k is the eigen- 
value of matrix (j^J corresponding to the energy of nth 
band for the wave vector k, with the eigenvector \n, k). 
Here Vi — \^jr* is the velocity of electrons (i — x,y,z). 
Taking the static limit of lj — ► 0, Eq. JJJ becomes 



n.k m^n 
2 



(8) 



where A„k = — i(nk\X7 k \nk) is the geometric vector po- 
tential and fnk is the Fermi-Dirac distribution function. 

The expression jHJ is similar to that obtained by Thou- 
less et alw& in the context of the Quantized Hall effect 
in two dimensions. In particular, when the band is com- 
pletely filled, the circulation of A n k or the flux of Berry 
curvature defined as tt n k — Vk x A„k over the first Bril- 
louin zone is equal to 2iri times an integeu22i2i called 
the Chern number v n . Expression © of the off-diagonal 
conductivity finally becomes 



fnk^nk = Y Y " n > 



(9) 



77 k 



where v n = ^ J f nii il nk d 2 k. As we shall see, this implies 
that the Hall conductivity is quantized when the Fermi 
level lies in a gap. 



The Chern numbers have the following properties: (i) 
the sum over all bands is equal to 0, (ii) if for some pecu- 
liar values of parameters of the model, the energy bands 
m and n cross each others at some point of the Brillouin 
zone, their respective Chern number obey the following 
conservation rule 



On 



V n )a 



(10) 



where the indices b, a refer to the values of the Chern 
number before and after the bands are touching each 
other We will check each of these conservation rules 
in the next sections. 

It should be noticed that the sum over k in Eq. 
runs over all occupied states. This may seem unusual as 
in metallic systems we expect that only the quasiparticles 
near the Fermi surface contribute to transport properties. 
However, it has been recently shown by Haldana2& that 
Eq. lO can be reconciled with this point of view and 
reduced to a sum over the states near the Fermi energy. 



III. NUMERICAL RESULTS 

This section is divided in two parts. In the first one, 
the band structure is calculated for different angles 9 and 
coupling constants Jo. The evidences for two different 
regimes of couplings are given, illustrated by the varia- 
tion of Chern numbers. In the second one, we focus on 
the transverse conductivity at zero temperature. It is 
shown that the chirality may be a relevant quantity for 
describing the transverse conductivity (a oc x) m some 
ranges of coupling only for a half-filled band. 



A. Band structure and associated Chern numbers 

The energy spectrum of the Kagome lattice in the ab- 
sence of exchange interaction, is characterized by one 
flat band at E = 2t and two dispersive bands, which 
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touch each other at the k-point K of the Brillouin zone 
as shown in Fig. |2 This peculiarity has already been 
addressed in the context of spin models as well as for the 
metallic systems2£. 

Once the coupling constant Jo is nonzero, the energy 
spectrum splits in two parts due to the spin dependant 
potential. For very large values of Jo, the spectrum is 
divided into two groups of three bands, which is the case 
studied by Ohgushi et alm&. In between, for nonzero 
but not too large values of Jo, the spectrum can only 
be computed numerically, except for the special cases of 
9 = or 7Tj or for general 9 at high symetry points. 
For these intermediate values of Jo , the splitting of the 
spectrum depends qualitatively on two mechanisms. 

First, the coupling Jo gradually separates each group 
of three bands taking them degenerate for Jo = to fully 
separated for infinite Jo. Second, within each group of 
three bands, point like degeneracies are lifted (see Fig. [21 
points T and K) when switching on J , and finally re- 
stored for Jo — > oo (see Ref . IT3l . The numerical calcu- 
lation of the spectrum shows that it is either gapless for 
small values of Jo, or has gaps for higher values. 

When gaps open, it always occurs at the M point of 
the Brillouin zone. This allows to compute analytically 
the critical value of the coupling as a function of chirality 
parameterized by the angle 9, 



Using Eq. ifTT)) we can distinguish between two different 
regimes (Fig.EJ depending on the value of Jo as compare 
to J c (9). These regimes are characterized by a particular 
organization of the band structure as shown in Figs. GJ 
(b) for Jq — t andEl(c) for Jo = 2t and different Chern 
numbers. For Jo > J c {9), the Chern numbers associated 
to each band are given by —1, 0, 1, 1, 0, —1 from the 
lower to the upper band. Those numbers were obtained 
by Ohgushi et alw^ in the limit of Jo — > oo. When Jo < 
J c , the Chern numbers associated to this configuration 
are —1, 3, —2, —2, 3, —1 given in the same order as in 
the previous case. The global sum rule J2n v n = is 
always verified as well as the local one around the critical 
coupling where two pairs of bands are crossing each other. 
It reads as 

v-i + ^3 = 1 and 1/4 + v 5 = 1 (12) 

below and above the critical coupling J c , respectively. 
The corresponding phase diagram is shown on figure 03 

B. Off- diagonal conductivity at T = 

The off-diagonal conductivity is computed using 
Eq. at T — for different values of the chirality 9 
and the amplitude of the exchange interaction J . As 
shown in Fig. Q] (the corresponding energy spectrum is 
presented in Figure 0(b) and 0(c)) at T = 0, the Hall 




FIG. 2: (a) Energy spectrum of 114t calculated for Jo = 0. 
Each band is twice degenerate due to the spin degeneracy, 
(b) Energy spectrum calculated for Jo — t and 9 — 7r/3. (c) 
Energy spectrum calculated for Jo = It and 9 = n/3. The 
critical value J c is equal to ktjy/l « 1.51 t 



conductivity is quantized when the Fermi level is lying 
in the gap. The value of the plateaus depends explic- 
itly on the Chern number of the filled bands. When the 
Fermi level is in a band, the sum can be reformu- 
lated in terms of a sum over the Fermi level as shown 
by Haldane^Sr recently. One has to note that when the 
Fermi level and 9 are fixed, any variation of Jo around 
J c (9) may induce a jump in the Hall conductivity. For 
instance with a filling factor p = 1/3 (Fig.QJa), the Fermi 
level always lies in a gap but the off-diagonal conductivity 
jumps from —e 2 /h for Jo > J c to 2e 2 /h for Jo < J c . 

A similar jump can be obtained of 9 is changed, while 
Jo is constant: such a change in 9 can for example 
be induced by application of an external magnetic field 
perpendicular to the kagome plane. In this case, large 
change of Hall conductivity (even sign change) can be 
induced by magnetic field. 

This is not the case for p — 1/2. For that filling fac- 
tor, the Hall conductivity varies smoothly when J passes 
through the point Jo = 3i/2 which does not depend on 
9 (Fig. GJb). When J < 3i/2, the Fermi level crosses 
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FIG. 3: The phase diagram in Jo — 9 variables. The black 
line is given by Eq. Ill IB and corresponds to the vanishing gap 
between the bands 2 and 3 (the gap between the bands 4 and 
5 is also vanishing). 
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the bands three and four resulting to a nonzero Hall con- 
ductivity. When J > 'it/2, the Fermi level is in the gap 
separating the band three and four and the off-diagonal 
conductivity is given by the sum off the Chern number 
of the first three bands which gives zero. 

Fig shows the variation of the Hall conductivity as 
a function of Jo for filling factors 1/4 (Fig, EJ-c) and 1/5 
(Fig- 13-d): in both cases the Fermi level is not in a gap 
and Hall conductivity does not exhibit any jump. 

From these comments, it is clear that the chirality is 
never a relevant parameter for the Hall conductivity when 
the Fermi level lies in a gap, because a xy is then quantized 
and obviously, not proportional to the chirality. 

Conversely, it is possible to choose a filling factor such 
that the Fermi level is not within a gap. In Fig. we 
present a xy versus 9 for the filling factor p = 1/2 (a) 
and (b) and different values of Jo. For this filling factor, 
the Fermi level is not in a gap if Jo < 3t/2. Scaling 
the conductivity by its maximal value in each case shows 
that it is proportional to chirality in a large range of 
coupling Jo , but also for large values of the chirality. This 
is somewhat unexpected as the previous studies relating 
v X y to x were restricted to small values of However, 
for different values of the band filling this is no longer true 
as shown on figureEJfor p — 1/4 (j^J-c) and p — 1/5 ®d): 
for these filling factors, a xy is sensitive to the variation of 
the Chern numbers when J varies. In fact we found that 
the Hall conductivity is proportional to chirality only for 
p = 1/2 and small enough Jo . 



FIG. 4: Off-diagonal conductivity evaluated at T = as a 
function of £/ for 9 — 7r/3 and (a) Jo = 2 t (b) Jo = t. In 
both cases, a xy is quantized when the Fermi level is lying in 
a gap. The values of the plateaus depend strongly on the 
topology of the energy spectrum. 



IV. DISCUSSION 

The model studied in this paper cannot be directly 
applied to materials, in which the Anomalous Hall Effect 
attributed to the spin chirality has been observed. 

However, some interesting features have been obtained 
in this model. This work shows that transverse con- 
ductivity may be proportional to chirality in a strong 
coupling regime, extending previous results obtained in 
the weak coupling case& but this occurs only for a half 
filled band. It is also shown that the dependence of 
a xy is not only determined by the chirality but also de- 
pends strongly on the exchange coupling parameter and 
on the band filling factor. a xy can present a large vari- 
ety of behaviors: non-monotonic variation, sign change 
or plateaus can be observed. Thus, even if the underlying 
magnetic phase chirality is the main origin of this intrin- 
sic anomalous Hall effect, it is far from being entirely 
determined by the chirality. 

We suggest that in real systems, variation of coupling J 
can be induced by temperature: the temperature acts on 
the system in many different ways. It acts on the mag- 
netic configuration of the localized spins Si by changing 
the amplitude of the magnetization, and on the position 
of the Fermi level in the band through the Fermi-Dirac 
distribution. 

We assume that the magnitude of each local moment 
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FIG. 5: Off-diagonal conductivity (J xy {J) evaluated at T = 
for a filling factor p — 1/3 (a), p = 1/2 (b), p = 1/4 (c) and 
p — 1/5 (d) and 6 = 7r/3 (continuous lines) and 6 = 7r/4 
(dashed lines). 



follows the mean field relation 



M{T) = | S» 1 = Mq\ 1 



T 



(13) 



where T c is the critical temperature. This can be de- 
scribed by introducing in eq. 1 a temperature dependent 
coupling constant 



J(T) = JoWl 



T 
TV 



(14) 



where Jo = JMq is the zero temperature exchange con- 
stant and Jo and T c are supposed to be independent. 

Consequently, decreasing the temperature from T c to 
T = is equivalent to increase the exchange coupling 
from J = to J = Jo- This increase of J can be consid- 
ered as done at zero temperature since in a large range 
of parameters, the additional effect of temperature (i.e. 
the T-dependance of the Fermi function in Eq. (8)), is 
significant only when the Fermi level is very close to a 
band edge. 




FIG. 6: a xy /ao with ao = m&x(\a xy \) versus 8 angle, calcu- 
lated for a filling factor p = 1/2 (a) and (b) and p = 1/4 (c) 
and p = 1/5 (d) and different values of Jo. The continuous 
line represents the chirality given by The Hall conduc- 
tivity is represented by the dashed lines for different values of 
the parameter Jo. We remark that a xy changes sign when Jo 
tends to 0.6i. That change is related to the position of the 
Fermi level within the band but not associated to change of 
the Chern numbers. It is remarkable that the conductivity 
is more or less proportional to the chirality in a large range 
of couplings, up to Jo w t as well as for large x values. For 
p = 1/4 and p = 1/5, one can show that the off-diagonal con- 
ductivity is not proportional to the chirality even for small 
values of Jo. 



Qualitatively, this means that when temperature is de- 
creased, one moves vertically from the bottom to the top 
of the phase diagram in Fig. if Jo is large enough, 
the temperature decrease produces at some peculiar tem- 
perature T* such that J(T*) = J C (T = 0), a change in 
conductivity due to the change of the Chern numbers. 
From this observation it follows that a xy can show abrupt 
changes with temperature. It is worth noting that dur- 
ing that process, in the frame of our model, the magnetic 
texture characterized by its 9 angle has not moved. Of 
course, a change in 8 can also be responsible for a sign 
change. 

To be more specific, we consider several examples pre- 
sented in Fig. [3 
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FIG. 7: Off-diagonal conductivity a xy (T) for a filling factor 
p = 1/2 (a), p = 1/3 (b), p = 1/4 (c) and p = 1/5 (d) and 
8 — 7r/3 (continuous lines) and = tv/4 (dashed lines). 



We start with the filling factor p = 1/2, the critical 
temperature T c = O.lt and J(T = 0) = 2t (Fig. da). 
With this choice of parameters, the Fermi level lies in 
a gap when Jo > 3 i/2, which means that at low tem- 
peratures (i.e., for large Jo), the system is insulating. 
When the temperature is increased, a xy remains equal 
to zero untill T reaches the temperature T\ for which 
J(Ti) = 3 i/2. Above this temperature, the gap closes, 
making the transverse conductivity increasing continu- 
ously. It reaches an extremum which is a function of 
the texture angle 6. Just below T c , a xy changes sign 
because of a subtle balance between the states giving 
positive and negative contributions to the conductivity, 



namely, the states from bands 3 and 4 (negative contribu- 
tion) and the states from band 2 (positive contribution). 
These contributions may be globally explained by the 
associated Chern numbers of these bands, at small Jo : 
^3, v a — —2 and = 3. 

For the filling factor p = 1/3 (Fig. db), the behav- 
ior is completely different. The transverse conductivity 
is quantized and finite at zero temperature as the Fermi 
level is located in a gap. When the temperature increases 
and reaches the value defined by J{T±) = J c (0), for 
which gaps are finite but sufficiently small to autho- 
rize the interband processes, a xy continuously increases. 
Upon increasing temperature, the conductivity increases 
and then changes its sign. This change of sign is not of 
the same origin as we considered for p = 1/2. In the first 
case, the Chern numbers are fixed but the temperature 
changes the balance between the weights of each band. 
In the second example, bands 3 and 4 get crossed thus 
changing their Chern numbers. The variation of Chern 
numbers produces the sign change in the transverse con- 
ductivity. Due to thermal fluctuations, the extremum of 
conductivity does not reach its quantized value of 2 (in 
units of e 2 /h) but saturates at 1.8. The conductivity fi- 
nally decreases down to zero when T reaches T c , as it 
should be. 

It is also possible to choose the filling factors that place 
the Fermi level at zero temperature within a band. In 
these cases, there is no generic behavior, and the angle of 
the spin texture plays an important role. Some results are 
shown on figured for the different filling factors p = 1/4 
(c) andp = 1/5 (d). 

Thus, these results show that in the mean field approx- 
imation a large variety of behaviors are possible in our 
model. The next step would be to study a model closer 
to the experimental situation of pyrochlores, where the 
variation of magnetic structure with temperature and ap- 
plied field has been studied by neutrons experiments^, 
allowing to take into account the real cristallographic and 
magnetic structures of these systems. 
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